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Abstract 

3C-like proteinase of severe acute respiratory syndrome (SARS) coronavirus has been demonstrated to be a key target for drug 
design against SARS. The interaction between SARS coronavirus 3C-like (3CL) proteinase and an octapeptide interface inhibitor 
was studied by affinity capillary electrophoresis (ACE). The binding constants were estimated by the change of migration time of the 
analytes in the buffer solution containing different concentrations of SARS 3CL proteinase. The results showed that SARS 3CL 
proteinase was able to complex with the octapeptide competitively, with binding constants of 2.44 x 10 4 M _1 at 20 °C and 
2.11 x 10 4 M _1 at 37 °C. In addition, the thermodynamic parameters deduced reveal that hydrophobic interaction might play major 
roles, along with electrostatic force, in the binding process. The ACE method used here could be developed to be an effective and 
simple way of applying large-scale drug screening and evaluation. 

© 2005 Elsevier Inc. All rights reserved. 
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In late 2002, severe acute respiratory syndrome 
(SARS) * 1 infected more than 1000 people, mostly in Chi¬ 
na, Hong Kong, Taiwan, and Canada, and caused deaths 
due to the infection by a novel human coronavirus [1-3]. 
Then the sequence of the complete genome of this novel 
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human coronavirus, SARS-associated coronavirus 
(SARS-CoV), was determined. It contained 29,727 nucle¬ 
otides in length and 11 open reading frames, with a gen¬ 
ome organization similar to that of other known 
coronaviruses, including two human coronaviruses, 
HCoV-OC43 and HCoV-229E. The phylogenetic analy¬ 
sis and sequence comparisons indicated that SARS-CoV 
was not closely related to any of the three previously 
known groups of characterized coronaviruses [4,5]. The 
SARS-CoV main proteinase (33.8 kDa), also called the 
3C-like (3CL) proteinase, was essential for the life cycle 
of SARS-CoV, playing a pivotal role in mediating viral 
replication and transcription functions through extensive 
proteolytic processing of two replicase polyproteins: ppla 
(486 kDa) and pplab (790 kDa) [6,7]. 

Because of the control of the activities of the corona¬ 
virus replication complex and the availability of 
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numerous inhibitors of other coronavirus proteinases 
and efficient expression, the 3CL proteinase was a pre¬ 
ferred target for the task of discovering drugs against 
SARS. Hilgenfeld’s group determined the crystal struc¬ 
tures of human coronavirus (strain 229E) main protein¬ 
ase and a substrate analogue hexapeptidyl chloromethyl 
ketone (CMK) inhibitor complex of porcine coronavirus 
(transmissible gastroenteritis virus, TGEV) main pro¬ 
teinase. The results revealed a remarkable degree of con¬ 
servation of the substrate-binding site, and a homology 
model for SARS CoV main proteinase was constructed 
[6]. The crystal structures of the SARS-CoV main pro¬ 
teinase at different pH values and binding to a specific 
inhibitor with a covalent-bonded substrate analogue 
were reported by Rao’s group, providing insights to 
the substrate binding site and a structural basis for ra¬ 
tional drug design [7]. The virtual screening tested by 
Xiong et al. [8] indicated that 73 available proteinase 
inhibitors in the MDDR database might dock into both 
the binding pockets of the TGEV main proteinase and 
the SARS-CoV 3CL proteinase. A series of cell-based 
assays with SARS virus and Vero E6 cells, including en¬ 
zyme-linked immunosorbent assay (ELISA), Western 
blot analysis, immunofluorescence, and flow cytometry 
assays, were applied and reported by Wu et al. [9] to 
screen existing drugs, natural products, and synthetic 
compounds to identify effective anti-SARS agents, with 
the results showing that approximately 50 compounds 
were active at 10 pM and that 7 compounds were active 
at 3 pM. Based on the atomic coordinates of the SARS- 
CoV main proteinase, research was conducted by 
Chou’s group on KZ7088 (a derivative of AG7088) 
and the AVLQSGFR octapeptide to the enzyme, which 
were interacting with the active site of the SARS enzyme 
through six hydrogen bonds, with KZ7088 having a 
clear definition of the binding pocket [10]. Kuo et al. 
[11] prepared a peptide with fluorescence quenching pair 
(Dabcyl and Edans) at both ends of a peptide substrate 
and used this fluorogenic peptide substrate to character¬ 
ize the SARS main proteinase and screen inhibitors. 

Capillary electrophoresis (CE) may offer several 
advantages to the study of the interactions between 
small ligands and biomacromolecules, including small 
sample consumption (nanogram), short analysis time, 
ease of automation, high efficiency, and high resolving 
power [12-14]. One company, Cetek, has developed 
and optimized the CE Assay, an advanced CE assay 
for high information binding [15,16]. The weak bindings 
(K b = 10 3 10 4 M -1 ) of macromolecules and small li¬ 
gands can be detected easily by affinity capillary electro¬ 
phoresis (ACE). Another advantage of the ACE method 
is that all of the interacting components can be studied 
in solution, thereby eliminating the possibility of dena- 
turation of conformational alteration of the biomacro¬ 
molecules such as protein and nucleic acid [12]. 
Various affinity interactions, such as drug-protein, pro¬ 


tein-protein, protein-nucleic acid, protein-carbohy¬ 
drate, peptide-antibiotic, enzyme-cofactors, lectin- 
sugar, antigen-antibody, and cyclodextrins (CDs)-enan- 
tiomer, have been investigated by ACE during recent 
years [17,18]. 

We demonstrated previously that the dimer of SARS 
3CL proteinase should be the biologically functional 
form and play a major role in catalysis [19], Inhibitors 
targeting the dimeric interface should be a new strategy 
to control the activity of SARS 3CL proteinase. Because 
one monomer of SARS 3CL proteinase binds to the 
other, specifically at the N-terminal interface [6,7], to 
compete with the dimeric interaction of the proteinase, 
an octapeptide interface inhibitor was designed accord¬ 
ing to the amino acid sequence of the N terminus of 
SARS 3CL proteinase. In the current study, a simple 
and effective CE method with very low sample consump¬ 
tion was developed and applied to study the interaction 
between SARS 3CL proteinase and the octapeptide 
interface inhibitor. The binding constants were mea¬ 
sured by the changes in the migration time of the octa¬ 
peptide inhibitor at different concentrations of SARS 
3CL proteinase in the running buffer solutions. The 
octapeptide was found to bind to SARS 3CL proteinase 
specifically and competitively. Thermodynamic analysis 
indicates that hydrophobic interactions play a major 
role, along with electrostatic force, in the binding 
processes. 

Materials and methods 

Reagents and materials 

SARS 3CL proteinase was expressed in Escherichia 
coli and purified according to published procedures [17]. 
The octapeptide (Ser-Gly-Phe-Arg-Lys-Met-Ala-Phe) 
was synthesized by solid-phase peptide synthesis using 
the standard 9-fluorenylmethoxycarbonyl/fer/-butyl 
strategy and was purified by reverse-phased HPLC. Bo¬ 
vine serum albumin (BSA, MW = 67,000) and ovalbumin 
(OVA, MW = 43,000) were purchased from Fluka. All 
other reagents were of analytical grade without further 
purification. Deionized water was used throughout. 

Apparatus 

CE was performed with the Beckman P/ACE 5000 
system (Fullerton, CA, USA). The uncoated fused silica 
capillary (Beckman, 57 cm x 50 pm i.d., 50 cm to the 
detector) was used. The capillary chamber temperature 
(±0.1 °C) was controlled by forced liquid cooling. A 
run voltage of 15 kV in the normal polarity mode was 
applied. UV detection was performed at 214 nm. The 
GS-15R multipurpose refrigerated centrifuge (Beckman) 
was applied to clear any trapped air bubbles in the 
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running buffer solutions containing SARS 3CL protein¬ 
ase before use. 

Solution and sample preparation 

The concentrations of SARS 3CL proteinase stock 
solution were determined from UV absorbance at 
280 nm (concentration/mg ml -1 = A 2 so x 1.04). The 
stock solutions of octapeptide, BSA, and OVA were pre¬ 
pared by accurately weighing defined amounts of the 
compounds, dissolving them in a defined volume of 
deionized water, and then storing at 4 °C. The running 
buffer in the CE experiments was 25 mM sodium phos¬ 
phate (pH 7.4). The investigated octapeptide and SARS 
3CL proteinase solutions in 25 mM sodium phosphate 
buffer (pH 7.4) were prepared by diluting the defined 
volume of the stock solutions to the constant volume 
with buffer solution. The constant concentration of 
octapeptide in the sample was 1.50 x 10 4 M. Dimethyl 
sulfoxide (DMSO) served as an electroosmotic flow 
(EOF) marker and was added to the sample at a concen¬ 
tration of 0.01% (v/v), which proved to have no influ¬ 
ence on SARS 3CL proteinase with the concentration 
less than 5% (v/v). The varying concentrations of SARS 
3CL proteinase in the running buffer were prepared 
from 0 to 10 pM. All solutions and deionized water were 
stored at 4 °C, being filtered through 0.22-pm cellulose 
acetate membrane filters (Shanghai Xingya Resin, 
Shanghai, China) and sonicated for 10 min to remove 
air from solution prior to use. SARS 3CL proteinase 
solutions were centrifuged for 5 min with 12,000 rpm 
at 4 °C to exclude air from solutions before use. 

Procedures for CE experiments 

New capillary was conditioned for 60 min with 0.2 M 
NaOH, for 30 min with water, and for 10 min with run¬ 
ning buffer. Between measurements, the capillary was 
flushed for 3 min with 0.2 M NaOH, for 2 min with water, 
and for 2 min with running buffer. The sample containing 
octapeptide and DMSO was introduced into the capillary 
by high-pressure injection (15 psi for 10 s). The electro¬ 
phoresis was carried out using phosphate buffer solutions 
containing different concentrations of SARS 3CL pro¬ 
teinase (0-10 pM). All separations were in triplicate. Rel¬ 
ative standard deviations (RSDs) of the migration times 
were calculated from a series of three experiments carried 
out with the same sample within 1 day. 

Results and discussion 

Determination of binding constants 

The binding constants (K b ) are the most fundamen¬ 
tally important parameters for biologically active mole¬ 


cules [12-14]. The principle of this method for the 
determination of binding constants is to exploit the 
changes in the migration time of the complexed sample 
in the background electrolyte containing the complexa- 
tion agent. After adding SARS 3CL proteinase to back¬ 
ground solution (BGS), the octapeptide inhibitor can 
interact with the proteinase. The equilibrium repre¬ 
sented by the following equations assumed that a 1:1 
complex is formed: 


OP + 3CLP ^ OP-3CLP 

(1) 

[OP-3CLP] 
b “ [OP] [3CLP] ’ 

(2) 


where OP, 3CLP, and OP-3CLP are the octapeptide, 
SARS 3CL proteinase, and complex, respectively. K b 
is the binding constant for Eq. (1). As pointed out by 
Yang et al. [20], the capacity factor k, which is used as 
an indicator to describe the relationship between reten¬ 
tion time and mass distribution equilibrium in chromat¬ 
ography, can be used here to describe the migration 
behaviors of the active compounds in CE. Moreover, 
the addition of SARS 3CL proteinase in the BGS is ap¬ 
plied as a pseudo-stationary phase. Based on our previ¬ 
ous work [21,22], the binding constant can be measured 
according to the following equation: 

k = ^^0+ K ^CLP}), (3) 

A'GmO 

where <p is phase ratio, [OP] s is the concentration of the 
octapeptide in stationary phase, and [OP] m0 is the initial 
concentration of the octapeptide without the addition of 
SARS 3CL proteinase in the mobile phase. 

As defined in chromatography, the capacity factor k 
is equal to (t - t 0 )/t 0 , where t and t 0 are the migration 
times of the sample and neutral marker, respectively. 
With the linear relationship of k versus the concentra¬ 
tion of SARS 3CL proteinase ([3CLP]), the binding con¬ 
stant is calculated from the slope (a) and intercept (b): 


The weak interactions between several anti-HIV com¬ 
pounds with HIV transactivation-responsive (TAR) 
RNA and BSA were successfully studied simultaneously 
and by quickly using this method with the binding con¬ 
stants of 10 3 to 10 4 M -1 [21,22]. 

To evaluate the pharmacological activity of inhibitor 
binding to SARS 3CL proteinase in a physiological sys¬ 
tem or simulated physiological system, the phosphate 
buffer system of 25 mM at pH 7.4 was used for the bind¬ 
ing study. The migration times of the octapeptide (t) and 
DMSO (/ 0 ) were measured in seven running buffer solu¬ 
tions containing different concentrations of SARS 3CL 
proteinase (0-10 pM). The experimental data and RSD 
values are summarized in Table 1. The reproducibility 
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Table 1 

Experimental data of the migration times and RSDs (n = 3) 


[3CLP] (x 1(T 6 M) 

20 °C 





37 °C 





t (min) 

RSD (%) 

to (min) 

RSD (%) 

k 

t (min) 

RSD (%) 

t 0 (min) 

RSD (%) 

k 

0 

6.573 

0.48 

8.050 

0.47 

-0.183 

4.704 

0.20 

5.898 

0.36 

-0.202 

0.953 

6.715 

0.10 

8.328 

0.10 

-0.194 

4.791 

0.22 

6.089 

0.12 

-0.213 

1.91 

6.821 

0.11 

8.511 

0.14 

-0.199 

4.838 

0.18 

6.180 

0.11 

-0.217 

3.81 

6.966 

0.13 

8.787 

0.22 

-0.207 

4.969 

0.10 

6.411 

0.26 

-0.225 

5.72 

7.137 

0.12 

9.108 

0.18 

-0.216 

5.083 

0.11 

6.603 

0.34 

-0.230 

7.62 

7.289 

0.17 

9.376 

0.23 

-0.223 

5.186 

0.10 

6.796 

0.14 

-0.237 

9.53 

7.407 

0.28 

9.600 

0.36 

-0.228 

5.305 

0.08 

7.008 

0.15 

-0.243 


Note. [3CLP], concentration of SARS 3CL proteinase; t, migration time of the octapeptide; t 0 , migration time of DMSO; k, capacity factor. 


of the migration times was very well based on the values 
of RSD less than 0.5%. A set of electropherograms of 
migration time shift of the octapeptide inhibitor versus 
SARS 3CL proteinase concentration change is shown 
in Fig. 1. A plot of k versus the concentration of SARS 
3CL proteinase at two temperatures based on the Table 
1 data is shown in Fig. 2. The binding constants of the 
inhibitor with the proteinase were determined four times 
at two temperatures and calculated according to Eq. (4), 
with the average values being 2.44 x 1 0 4 VI 1 
(RSD = 6.0%, n = 4) at 20 °C and 2.11 x 10 4 N'T 1 
(RSD = 15%, n = 4) at 37 °C. The RSD value at 37 °C 
was larger than that at 20 °C due to the instability of 
the higher capillary temperature. 
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Fig. 1. Electropherogram of the migration time shift of the octapep¬ 
tide versus SARS 3CL proteinase concentration change at 20 °C. 
Conditions: running buffer, 25 mM phosphate solution (pH 7.4); 
applied voltage, +15 kV; pressure injection, 10 s at 15 psi; wavelength, 
214 nm. Peaks 1 and 2 represent octapeptide and DMSO, respectively, 
and A, B, C, and D represent SARS 3CL proteinase concentrations of 
0, 0.953, 1.91, and 3.81 pM, respectively. The constant concentrations 
of octapeptide and DMSO in the sample were 1.50xl0~ 4 M and 
0.01% (v/v), respectively. 


We demonstrated previously that the dimer of SARS 
3CL proteinase should be the biologically functional 
form and plays a major role in catalysis. In addition, 
one monomer of SARS 3CL proteinase binds to the 
other one specifically at the N-terminal interface 
[19,23]. To compete with the dimeric interaction of the 
proteinase, the octapeptide inhibitor was designed 
according to the amino acid sequence of the N terminus 
of SARS 3CL proteinase. Because the concentrations of 
SARS 3CL proteinase were less than 0.2 mg mV 1 in the 
current experimental conditions, the main form of the 
proteinase was thought to be monomer. The binding 
constant of the octapeptide with SARS 3CL proteinase 
was measured to be 2.44 x 10 4 M _1 at 20 °C, and the 
dissociation constant of the dimer of the proteinase 
was estimated to be 100 pM [19]; thus, the octapeptide 
inhibitor can bind to the monomer of the proteinase 
competitively and can prevent the dimerization effec¬ 
tively. Moreover, the interactions between the octapep¬ 
tide and the two common proteins, BSA and OVA, 
were studied as a negative control at the same conditions 
to avoid the influence of nonspecific absorption. Com¬ 
pared with Fig. 2, the plots of BSA were similar 
(Fig. 3 A), with the binding constants being calculated 
as 8.01 x 10 3 M” 1 at 20 °C and 8.20xl0 3 M~' at 
37 °C, whereas the plots of OVA were random 
(Fig. 3B). The results showed that the octapeptide could 
interact with BSA weakly but could not bind to OVA 
and that there should be a specific binding between the 
octapeptide inhibitor and SARS 3CL proteinase. 

Thermodynamic studies of the interactions 

Small molecules bind to macromolecules with four 
types of interactions: H-bond, van der Waals, electro¬ 
static, and hydrophobic interactions. The thermody¬ 
namic parameters, enthalpy change (AH) and entropy 
change (ASj of the binding process, are important for 
confirming binding mode [24,25]. If AH does not vary 
significantly over the temperature range studied, the 
enthalpic contribution to the Gibbs free energy (AGj 
can be determined according to the following equation: 
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Fig. 2. Plot of k versus the concentration of SARS 3CL proteinase at 
two temperatures. 




Fig. 3. Plots of k versus the concentrations of BSA and OVA. (A) 
Interaction between octapeptide and BSA at 20 and 37 °C. (B) 
Interaction between octapeptide and OVA at 20 and 37 °C. 


In 



AH /1 1 \ 

~R\T[~Y 2 )^ 


(5) 


where K t and K 2 are the binding constants at tempera¬ 
tures T x and T 2 , respectively, and R is the gas constant 


Table 2 

Thermodynamic parameters of interactions between SARS 3CL 
proteinase and the octapeptide 


T (K) 

K b 

AH 

AG 

AS 


(x 10 4 M 1 ) 

(kJ moR 1 ) 

(kJ moR 1 ) 

(JmoR 1 KR'l 

293 

2.44 

-6.71 

-24.6 

61.1 

310 

2.11 

-6.71 

-25.7 

61.3 


Note. T, absolute temperature; K b , binding constant; AH, enthalpy 
change; AG, Gibbs free energy; AS, entropy change. 


(8.314 J mol 1 K 1 ). The free energy change is estimated 
from the following relationship: 

AG= -RTlrtK, (6) 


where K is the binding constant at the corresponding 
temperature. The entropy change can be determined 
from the following equation: 


AS = 


AH - AG 
T 


(V) 


Based on the values of binding constants, the thermody¬ 
namic parameters AH, AG, and A .S' were measured 
according to Eqs. (5)—(7), respectively. These values 
are summarized in Table 2. 

As shown in Table 2, AH values were negative un¬ 
der the experimental conditions, demonstrating that 
the binding reaction was an exothermic process. The 
sign and magnitude of the thermodynamic parameters 
associated with various kinds of interaction were char¬ 
acterized [24,25]. From the point of view of water 
structure, positive entropy was frequently taken as 
the evidence of hydrophobic interaction, and it was 
also shown that positive entropy and slightly negative 
enthalpy might be a manifestation of electrostatic 
interactions between ionic species in aqueous solution 
[25]. Based on the experimental data, we conclude that 
the hydrophobic interaction may play a major role, 
whereas electrostatic forces also contribute to the 
binding process of SARS 3CL proteinase and the 
octapeptide inhibitor, consistent with our previous 
work. 


Conclusion 

SARS therapy requires the development of antiviral 
compounds that effectively prevent or treat this disease. 
A simple and reliable CE method with very low sample 
consumption was applied to study the interaction be¬ 
tween SARS 3CL proteinase and an octapeptide inter¬ 
face inhibitor. The binding constant was determined 
by the changes in the migration time of the octapeptide 
inhibitor at different concentrations of SARS 3CL pro¬ 
teinase in the running buffers. The octapeptide was 
shown to bind to SARS 3CL proteinase competitively, 
with binding constant K b values of 2.44 x 10 4 M -1 at 
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20 °C and 2.11 x 10 4 M -1 at 37 °C. Analysis of the ther¬ 
modynamic parameters, enthalpy change (AH) and en¬ 
tropy change (AS), indicated that hydrophobic 
interaction might play a major role, along with electro¬ 
static force, in the binding processes. The ACE method 
should be widely applicable for large-scale ligand screen¬ 
ing against protein targets. 
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